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Using a two-dimensional linear least-squares method applied to Doppler 
radar data, we test the viability of determining F-factor remotely. The ultimate 
application of such an algorithm will be supplying real-time F-factor maps, derived 
from ground-based Doppler radars to air traffic control personnel and pilots. Data 
from NASA deployments to the MIT/Lincoln Lab TDWR testbed radar in Orlando 
in 1991 and 1992 along with NASA deployments to the NCAR TDWR testbed 
radar in Denver are examined. Preliminary analyses show that the two- 
dimensional method correlates reasonably well with in situ measurements. 

Several effects, independent of the method used, act to reduce the correlation to 
less than 1. These include time differences between radar and aircraft data, 
vertical misalignment between the aircraft and the radar beam, different spatial 
resolution scales between aircraft and radar data, inhomogeneous radar beam 
filling, noise in radar data that eludes filtering and phase lag between time and 
space due to low pass filtering of the aircraft data. In the final assessment, it 
appears that a shear-based F-factor algorithm is preferable to the currently 
implemented TDWR algorithms which lack any local shear es tima tes 
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Summer ’92: NASA flies instrumented 737 
through Denver and Orlando microbursts 
close to radar beam height. Compare in 
situ aircraft data to shear estimates using 
TDWR testbed radar. We estimate F by 


F=-U r GS- 

9 


w 

GS 


where U r is radial shear. U r is estimated 
by a two dimensional least squares fit of 
the radial doppler velocity. We estimate w 

by 


w=- 2U r Z 

where Z is the height of the beam. We are 
assuming local symmetry of the 
microburst. This estimate for w is used at 
all points where we detect divergence 
(U r >0). 
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Using weighted least squares, solve 


and get 


EL = 


U e = 


ArY,i 2 Wij 
SA Si jUijWjj 

£(A Sij)*W itj 


where 


U{j = radial velocity at i, j, 


A Sij = ri<f>j and 
Wij = weight. 
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We estimate the horizontal F factor by 


, = 1 dV 
h gdt 


»-• U-GS 

g 


and the vertical F factor by 

F __ w., 2U z Z 

v GS GS 

when we detect divergence and 



GS 


when we detect convergence. 
Total F is approximated by 


F=F h +F v . 
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Factors Affecting Individual Path Correlations 


Time differences between aircraft penetration and radar data 
collection. 

Vertical mismatch between aircraft flight path and radar 
beam; aircraft either above or below the beam. 

Aircraft distance from beam (too close or too far from radar). 

Spatial resolution scale differences; the aircraft resolves 
smaller scales than the radar. 

Inhomogeneous radar beam filling; scatterers may not be 
evenly distributed throughout the beam illumination volume. 

Small scale noise within radar data that eludes median 
filtering. 

Phase lag between time and space of aircraft data due to low- 
pass filtering. 
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b462 Flight Path 


TIME =92 07 24 01 27 16 



256 




25 ' 



Elevation (m) 


Radar 


Max Radar 


Elev Stats <b462 




259 






- 30.0 


Velocity Stats<b462 



260 




F„F Actor 


F_Fac tor Stats «b462 85th pet 'froth 


C\J 



Radius (km) 


261 




262 




Vert I cdl.F 

-0.05 -0.03 0.00 0.03 0.05 


Vertical_F Stats ib462 85th pet Truth 



263 





Estimated F Factor 

jp. 200-0. 170-0. 140-0. 11 0-0. 080-0. 050-0. 020 0.010 0.040 0.070 0.100 0.130 0.160 0.190 


F factor Scatter 



264 




?n 1 ' ' ' 

10.0 15.0 


Radius (km) 


265 






f _hor i z Stats «b462 85 th pet fruth 


(N 



Radius (km) 


266 



Vert lcdl_F 

-0.05 -0.03 0,00 0.03 0.05 


85th pet Tro 


Verticdl_F Stats*b462 



V> 

/ 

J 

T^V 

1 

l 

L 

1 1 L. 

'i 

1 1 


10.0 15.0 

Radius (km) 


267 



. 200-0. 1 70-0, 140-0. 11 0-0. OBO-0. 050-0. 020 0.010 0.040 0.070 0.100 0.130 0.160 0.190 


F factor Scatter 



268 


R&d i d 1 Velocity 


TIME =92 08 17 22 37 33 



269 


b555 Flight Path 


TIME =32 08 17 22 39 09 



270 



b555 Radial V 


TIME =92 08 17 22 39 09 



271 



b555 Radar Data 
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Additional Points 


Clean up the current data set; trim off regions where the 
aircraft path is not representative of radar beam path. 

Threshold (categorize) data and compare only those cases 
where either aircraft of radar-derived F exceeds 0.085 (or any 
other threshold). 

Examine F comparisons only within shapes. 

Current implementation does not constitute an operational 
algorithm; various problems peculiar to Doppler radars must 
be addressed for operational implementation. 

A method of translating F factor to hazard severity needs to be 
created. 
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Conclusions 


Path-specific statistics show that a least squares-based method 
faithfully reproduces significant shears experienced by 
penetrating aircraft. 

Microburst wind shears display local linear symmetry. 

A shear-based microburst detection algorithm identifies only 
areas containing a quantifiable performance loss. 

F contains potentially more information because it uses shear; 
the currently implemented algorithms discard all local shear 
estimates. 

We believe that a shear-based F-factor algorithm is preferable 
to the present TDWR headwind loss-based algorithm (These findings 
are similar to that being proposed for ITWS). 
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